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Abstract 

We present a novel class of real symmetric matrices in arbitrary dimension d, linearly dependent 
on a parameter x. The matrix elements satisfy a set of nontrivial constraints that arise from asking 
for commutation of pairs of such matrices for all x, and an intuitive sufficiency condition for the 
solvability of certain linear equations that arise therefrom. This class of matrices generically violate 
the Wigner von Neumann non crossing rule, and is argued to be intimately connected with finite 
dimensional Hamiltonians of quantum integrable systems. PACS 71.10 Fd, 2.30 Ik, 2.10 Yn 
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Introduction We present a novel class of real symmetric matrices in arbitrary finite 
dimensions d. These matrices are linearly dependent on a parameter x, which plays the role 
of an interaction constant for a quantum system. The matrix elements satisfy a set of linear 
as well as non linear constraints that are derived. Each matrix a{x) ( as in Eq(0)) of this 
class has commuting partner matrices j3{x) ( as in EqQ), also linearly dependent on x, and 
having d + 1 independently assignable real parameters. This class of matrices generically 
exhibits level crossings as x is varied, and is intimately connected with finite dimensional 
Hamiltonians of integrable systems with dynamical symmetries. 

The study of linearly parameter dependent matrices is very popular in the context of 
quantum chaos [H 0, y, where it models the change in universalty class of quantum 
systems, or the temporal variation of correlations within the same class via the matrix 
Brownian motion model of Dyson|^,0,0]. This work has a different goal, that of formulating 
constraints so that the system remains 'regular" despite mixing with another matrix. In fact 
such linearly parameter dependent matrices are prototypes of quantum integrable models. 
The well known case of the anisotropic Heisenberg model solved by Bethe's Ansatz Q 
possesses "infinite" higher conservation laws, every one of them being linear in the anisotropy 
parameter^. 

The immediate motivation for our work comes from a study of the parameter depen- 
dence of eigenvalues of blocks of finite dimensional matrices arising from cornpletely inte- 
grable models of interacting quantum systems, such as the Heisenberg model [q] or the 1-d 
Hubbard modeljlO]. The famous Wigner- von Neumann fWvN) non crossing rule for 
parametric evolution of eigenvalues in quantum mechanics is a fundamental guiding principle 
for understanding level crossings, with violations or exceptions being termed as "acciden- 
tal degeneracies". The term accident is used since usually there is no specific "space time 
symmetry" responsible for such a degeneracy. There is a general belief that dynamical sym- 
metries, i.e. operators dependent on the interaction parameter also lead to degeneracies and 
violations of the WvN rules, despite the lack of a general group theoretic argument of the 
type that space time symmetries allow|l2]. 

Dynamical symmetries occur in most integrable systems. Our interest is also in sharpen- 
ing the notion of complete or exact integrability in the context of finite dimensional systems. 
In classical mechanics we have a very clear statement about the meaning of complete inte- 
grability, namely that the number of degrees of freedom equals the number of conservation 
laws that are mutually consistent. The meaning of the term "degrees of freedom" is quite 
unambiguous. For example in the trivial case of a set of harmonically coupled oscillators, it 
is the number of oscillators so coupled. In quantum mechanics, we do have similar unam- 
biguous models, such as the Calogero-Sutherland-Moser svstem|l3l| in the continuum, or the 
Toda 1J| lattice, where there exists a natural definition of a degree of freedom, in complete 
parallel to the classical situation, essentially the number of "particle" type variables. How- 
ever, in the case of other quantum integrable models that arise in condensed matter physics, 
such as the Heisenberg spin chain or the 1-d Hubbard model, the situation is ambiguous: 
the number of spin flips, or particles is variable. These models are defined on a discrete 
lattice, and have a state space that is in general finite, and lead to finite dimensional matri- 
ces depending upon a parameter, say the spin space anisotropy or the interaction strength. 
Presented with a finite dimensional matrix arising by restricting such a model to a finite 
lattice, it is challenging to distinguish it from other matrices of the same dimension. With- 
out referring back to the defining parent models, it is generally impossible to recognize their 
being "integrable" , whatever that word implies! Indeed an extreme and skeptical view would 
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challenge the very notion of complete integrability in a finite dimensional setting. One may 
argue that there are always d commuting independent matrices for any given matrix H; one 
merely diagonalizes the matrix and in its eigenbasis, constructs the d projection operators 
P, = |j X j|, so that [H, P,] = = [P„ P,]. 

Our viewpoint in this work is that the parameter dependence of eigenvalues contains 
the essence of quantum integrability for such finite dimensional models. These lead to a 
violation of the WvN rules, and hence to Poisson statistics of the energy level separation. 
While the relationship between dynamical symmetries and levels crossings is not yet precisely 
established, there are several studies that indicate a deep relation ship between them. For 
example, the beautiful numerical work of Heilmann and Lieb in 1971 [ISj on the 1-dimensional 
Hubbard model on a six site ring, shows that after all the known space time and internal 
space symmetries are carefully extracted, the finite dimensional blocks of matrices labelled 
by the appropriate quantum numbers, displaying scant regard for the WvN rules, have a 
large number of level crossings as the interaction constant is varied. More recent work of 
Yuzbashyan et alfl^ has examined the detailed connection between these level crossings 
and the dynamical symmetry operators of the Hubbard model found in 1986 by Shastry[l7[, 
and provided considerable insight into this phenomenon. In particular, there is an explicit 
algebraic demonstration for d = 3 that dynamical symmetries definitely imply level crossings. 

Another related and prominent manifestation of integrability is that the energy level 
statistics of such models display Poisson type behavior that allows spacings to be arbitrarily 
small ^J. This is in sharp contrast with the level spacings of generic (i.e. nonintegrable) 
models that exhibit level repulsion as expected from the WvN rule, and follow one of the 
three typical behaviors relevant to their class of quantum systems- namely the Gaussian 
orthogonal, unitary or symplectic classes 3 • 



It remains however, to state explicit conditions on matrices in arbitrary finite dimensions, 
that could identify the proclivity for level crossings, and hence presumably define completely 
integrable cases without reference to a parent quantum model. This goal is achieved in this 
work, we present a set of conditions, and a class of matrices satisfying them in any dimension 
d. We find through examples that this class of matrices automatically leads to Poisson type 
statistics for the energy level separation, and also an abundance of level crossings. 

Our main results follow from asking for the conditions for two parameter dependent 
matrices to commute with each other for all values of the parameter. Upon imposing a very 
intuitive sufficiency condition of autonomy (Type I matrices as explained below), it leads to 
a set of constraints for each of the matrices. The partner matrix in the commutating pair is 
automatically also a member of the same class. 

The matrix equations: We consider real symmetric matrices in d dimensions, linearly 
dependent on a real parameter x through 

a{x) = a + xA (1) 

where a is a diagonal matrix with diagonal entries {oi, 02, aa} and A is a real symmetric 
matrixji^. In the matrix a we have d{d— l)/2 off diagonal variables Aij, d variables aj and 
a further d variables Aj = Ajj, in addition to the real variable x. In an identical fashion we 
consider another matrix 

P{x) = b + xB, (2) 

where the diagonal matrix b has entries {61, 62? ^cz} and B is another real symmetric matrix. 
Clearly at x = the matrices a and (3 commute. We now ask the question, what are the 
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conditions under which these commute for arbitrary values of x? This clearly leads to two 
independent conditions 

[a, B] = [b, A] and [A, B] = 0. (3) 

The first set of d{d— l)/2 conditions, are expressible in terms of an antisymmetric matrix 
Sij as follows: 

S^, = ^^ = J^ for ^<J. (4) 
ai — aj bi — Oj 

The second set of equations can be written compactly using the Eq(^ in terms of 
where Pij = 4^3^- The d{d — l)/2 equations® can be written as 

r,,H = F,,[/5], (6) 

where is obtained from the same formula Eq(0), with (a^, Aj) ^ {bj, Bj). Although it 

is possible to symmetrize Yij in i&j, for future use it is better to use the present unsymmetric 
form. As expected the Equation® is symmetric in a and f3. At this point, we disturb this 
symmetry, and rearrange terms so that we can solve for /3 if a were given. This implies with 

{b A a)ij^i = [{biaj - bjat) + {bjai - kaj) + {kai - biai)] 

6, = - B, = p,,{h - b,) + -4- E %^(& A a),,,. (7) 

This set of equations is linear in the 2d variables {&i,&2; --bd} and {Bi, B2, B^, ...,Bii}, and 
provides conditions for the matrix (3{x) to commute with a{x) for all x. The offdiagonals 
of the matrix B namely Bij are already fixed by Eq® once the b'jS are picked in terms of 
the Sij. For a given a{x) then, we have the freedom of choosing these 2d variables subject 
to the d{d — l)/2 constraints Eq®, which for large d, are many more than the number of 
variables available. These constraints are best expressed in terms of a "triangle law" for any 
three distinct indices k: 

C{i,j,k) = ^i,, + ^j,k + ^k,i = 0. (8) 

Note that C is antisymmetric under exchange of any pair of arguments. Since these con- 
straints are linear in b'jS we collect the coefficients and write 

C{i,j,k) = iJ,{i;j,k)bi + iJ,{j;k,i)bj + iJ,{k;i,j)bk+ J2 ^{hhhk)bi. (9) 
/i is defined below in Equation (fT?l|) and 

^ii;hj,k) = ^^ + ^^ + ^^. (10) 

^jik ^k,i 

The solutions of Equations® may be classified as being of two types. Type I solutions 
correspond to requiring the coefficients of every br in Equation® to vanish individually. 
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Such a choice is sufficient without being necessary, and gives us constraints on the S matrix 
all by itself. It also leads to autonomous constraints on a, i.e. constraints involving the 
variables aj, Aj, Sij, but not the j3 variables. Type II solutions include all other possibilities, 
where the coefficients of br are not all vanishing. These include the trivial solution a = f5, 
and are less interesting as such. Type I solutions give us families of matrices a as we show 
below, that are definable in terms of relations involving S"s themselves, and between these 
and the remaining variables. 

Continuing our study of Type I solutions, we therefore equate the coefficients of b's 
individually to zero, giving for all "^Cs distinct triples i,j, k the three index formulas, 

fi{i]j,k) = plus cyclic permutations of fc. (11) 

Likewise we get '^€4^ four index formulas for each distinct quadruple of indices, z/(/; i, j, k) = 0, 
or rearranging a bit: 

Si,iSijSj^kSk,i + Sj^iSi^kSk,iSij + Sk^iSi^iSijSj^k = 0. (12) 

Note that the expression on the LHS is fully symmetric in the four variables. Although this 
relation involves quartics in the S"s, inspection shows that for the case of real symmetric 
matrices, a considerable simplification occurs and this constraint can be written in terms 
of bilinears if one inverts the matrix elements ( not the matrix itself!) of S and defines 

= (13) 

The vanishing of u in Equation(|12p can be restated as the ( fermionic Wick's theorem type) 
requirement that the totally antisymmetric symbol (p vanishes for all distinct quadruples 
i, J, k, I 

4'ii,j,k,i) = Ri,jRk,i — Ri,kRj,i + Ri,iRj,k = 0. (14) 

For large d, the four index constraints Eq p4|) are ~ d'^/ 4! in number, representing a huge 
overdetermination since the number of variables available, namely the Sij are only ~ d'^/2 
in number 2l|. Fortunately these identities are not all independent, and there exists an 
important extra identity relating a set of five indices that can be proved. For any distinct 
set of five indices k, /, m, we can easily see that 

4'{i,j,k,l)Rl,m = —4'{i,j,l,m)Rl,k + 4'{i,k,l,m)Rl,j — (j) {j ,k ,1 ,m) Rl ,i ■ (15) 

Using this five index identity we can show that the number of independent quadruple rela- 
tions are only {d — 2){d — 3)/2 in number. These may be chosen to be 



(1,2,3,4) 



0, 0(1,2,3,5) = 0, ... 0(l,2,d-l,d) = 0. (16) 



From these equations, we can satisfy all others using Eq (fT3j) repeatedly. As an example 
consider 0(3,4,5,6), we merely multiply by R^^i so that on using Eq lfTHjl . we find 

0(3,4,5,6)-R3,l = 0(3,1,5,6)-R3,4 - 0(3,1,4,6) -^3,5 + 0(3,1,4,5)-R3,6, (17) 

and each of these can be similarly processed further, e.g. 

0(1,3,4,6)-Rl,2 = 0(1,2,4,6)-Rl,3 " 0(1,2,3,6) -^1,4 + 0(1,2,3,4)-Rl,6 (18) 
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and reduced to forms involving 0(i,2.i,m) which vanish according to our hst of Eq()16|). 
We now turn to the study of fi{i',j, k) which is written compactly as 

M^;j,A:) = F,,-r,,-%^ (19) 



involving the same variable Yi j that we encountered in Equation© j22]. The function ix 
is antisymmetric in j, k as it stands. It can be antisymmetrized in all three variables by 
defining k) = (a^ — aj){ai — ak)fi{i;j, k), which satisfies 

k) = (a, - aj}[Ak '- — + [aj - ak}[Ai ^ — 

+ (afc - a,)[Aj - ^1^1^] - (ai - aj){ai - a^) J2 ^i.r ( ^ - ^) • (20) 

k,i \ i,j i,k / 

We comment on some important properties of these formulae. The first part of Eq ()2()|l 
consisting of three terms is explicitly antisymmetric in the three indices. The last line is 
antisymmetric in j, k but not manifestly so in i, j. We add to it a term with z, j exchanged, 
leading to 

S' S S' S S' s 
fx{i,j,k) + fx{j,i,k) = (ai - aj){ak - ar)[ ''"^ '''^ + ^'^ ''^ (21) 

rjti,j,k ^j^f' 

but it vanishes on using the four index formula Eq (jl2|l . whereby fi{i,j,k) is antisymmetric 
in all three indices. 

One more beautiful property of is: 

/i(z; J, k) - /i(z; J, I) = Y,, - F,, + %^ + = /i(z; I, k), (22) 

where the last line follows again from the use of the four index formula Eq (fT^ . Thus the 
total number of independent constraints of the /i type are only (d-2) in number (2 lii]. instead 
of the apparently huge number '^C^- We may choose them most simply as 

/i(l; 2, 3) = 0, .../i(l; 2, - 1) = 0, /i(l; 2, d) = 0. (23) 

The reader can verify that all other constraints of the type fi{i,j, k) = are obtainable from 
these d — 2 by using Eq ()22j) and the antisymmetry of fi{i,j, k). 

Constraints, variables and consistency: We now recount the number of available 
variables versus the constraints and show how generic matrices of Type I can be constructed. 
Firstly we construct the S matrix satisfying the {d — 2){d — 3)/2 constraints Eq()16p: this 
can be done by assigning arbitrary values to {2d — 3) matrix elements of 5" and computing 
the rest from the constraints. A particularly convenient choice is to assign values to Sij 
and S'2j for {2, 3} < j < whereby the constraints Eq(fTH|) reduce to linear equations for 
the other matrix elements. 

Having determined the S matrix, we next construct the a matrices. A straightforward 
strategy is to assign arbitrary values to the d + 2 variables {%}, ^i, ^2, and then using the 
d — 2 constraints Eq lj^Hj) as linear equations for the remaining Aj. Thus the total number 
of freely assignable variables for constructing the a matrix is 3(i — 1. 
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We can next determine a f3 matrix by assigning arbitrary values to the set of d + 1 
variables {bj} and one of the B/s (say Bi), and using the linear Equations© for ^j^i to 
determine the rest. Thus c? + 1 is the number of independent (3 type matrices of Type I for 
a given a. 

By construction the resulting f3 and a matrices commute for all x. It is further clear that 
the Equations® are satisfied identically. Hence it follows that fi{i;j,k) = for all triples 
is satisfied whether we use Yij{a) or Yij{f3). This guarantees that starting with a matrix 
a{x) of Type I, the resulting matrix i3{x) is automatically of Type I. 

We noticed that in addition to the solutions presented above, there are some beautiful 
special solutions of these constraints Equation p4|) belonging to the class of Toeplitz matrices. 
From inspection and using various addition theorems, it is clear that the following class of 
R's satisfy this constraint identically: 

Rij = - or Rij = sm{K{i) - ^(j)), (24) 

where K{j) is an arbitrary function of its argument. 

We plan to return to the problem of Hermitean as well as symplectic matrices in a future 
workj23|. One expects similar results to the ones presented here, but with more elaborate 
constraints. The results in ci = 3 are contained in our present ones, and correspond to 
the simple case of requiring /i(l; 2, 3) = 0; this result is sufficient to make the discriminant 
vanish at points in x. For higher d the analysis of the discriminant is more difficult, however 
examination of several examples of Type I matrices leads us to conjecture that these always 
lead to level crossings\2^. A direct algebraic proof involves examining the condition for 
vanishing of the discriminant of the matrices, and is currently being pursued. 

An important issue concerns translation invariance of the results for Type I matrices. It 
is not obvious that a matrix of Type I remains so if we shift x by a constant and absorb 
the change into the diagonal piece a by rediagonalizing using a suitable orthogonal trans- 
formation. We have verified that the results do possess this translation invariance in the 
parameter x, both numerically ( for small d) and analytically. A detailed proof using the 
parameter derivatives of the constraints Eqs ()23ll2j) . using the Pechukas flow equations |2^ 
is essentially complete, and will be published separately. In future work we hope to address 
several physical models using these constraints ( and some obvious variants), to check for 
their compliance. 

I acknowledge support from the NSF grant DMR 0408247. 
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